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Abstract 

This paper deals with the question of conditional sampling and prediction 
for the class of stationary max-stable processes which allow for a mixed 
moving maxima representation. We develop an exact procedure for conditional 
sampling using the Poisson point process structure of such processes. For 
explicit calculations we restrict ourselves to the one-dimensional case and use 
a finite number of shape functions satisfying some regularity conditions. For 
more general shape functions approximation techniques are presented. Our 
algorithm is applied to the Gaussian extreme value process and the Brown- 
Resnick process. Finally, we compare our computational results to other 
approaches. 
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1. Introduction 

Over the last decades, several models for max-stable processes have been developed 
and applied, including times series models and stochastically continuous processes. In 
view of the wide range of potential applications of max-stable processes for modelling 
extreme events, the question of prediction and conditional sampling arises. Davis and 
Resnick (|S|, proposed prediction procedures for time series which basically aim to 
minimize a suitable distance between observation and prediction. Further approaches 
for time series or random fields have been rare for a long time, apart from a few 
exceptions. Cooley et al. [3] introduced an approximation of the conditional density. 
Recently, Wang and Stoev 22 proposed an exact and efficient algorithm for conditional 
sampling for max-linear models 

Zi = max dijYj, i = 1, . . . , n, 

where Yj are independent Frechet random variables. Dombry et al. [S] presented an 
algorithm for conditional simulation of Brown- Resnick processes based on more general 
results on conditional distributions of max-stable processes given in [7]. 
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Here, we will consider stationary max-stable processes with standard Frechet mar- 
gins that allow for a mixed moving maxima representation (see, for instance, |17| . |20j). 
Let {^l,J-,P) be a probability space and F : {^,J-) — > {G,G) be a random function 
such that E(/jj^ F(x) dx) = 1. We assume that G is a countable set of measurable 
functions f : R'^ ^ [0,oo) and G = 2^. 

Then, we consider the stationary max-stable process 

Z{t)^ max uf{t-s), teW^, (1.1) 
where 11 is a Poisson point process on S" = M'^ x (0, oo) x G with intensity measure 
A{AxBxC) = \A\ ■ Pf{C) • / 4> A e B"^, B e B n {0,oo), C e G, (1.2) 

Jb " 

and Pf is the probability measure belonging to F. 

We aim to sample from the conditional distribution of Z given Z(ii ),..., Z(t„) 
for fixed ii, . . . , t„ e M"^. As Z is entirely determined by the Poisson point process 
n, we analyse the distribution of 11 given some values of Z. The idea to use a 
Poisson point process structure for calculating conditional distributions has already 
been implemented in the case of a bivariate min-stable random vector |23| . A very 
general Poisson point process approach was recently used by [7] yielding formulae for 
conditional distributions in terms of the exponent measure. Some of these results are 
independently found here, as well. 

The paper is organized as follows. In Section [2l we introduce a random partition 
of n into three measurable point processes. This partition allows to focus on the 
critical points of 11 which determine Z{ti), . . . , Z{tn). Similarly to [22_ and [7^, we will 
call realisations of these point configuration scenarios and figure out the conditional 
distribution of these scenarios coping with the problem that we work on events of 
probability zero (Section[3]). In Section|4l the conditional distribution of 11 is calculated 
explicitly for the case d — 1 and some regularity assumptions on a finite number 
of random shape functions. In Section [5l the results are applied to the Gaussian 
extreme value process 119] and compared to other algorithms. Section [6] deals with 
an approximation procedure in the case of a countable and uncountable number of 
random shape functions. A prominent example, the Brown- Resnick process [T], is the 
matter of a comparison study for different algorithms in Section [T) In the last section, 
we give a brief overview of the results for a discrete mixed moving maxima process 
restricted to pZ'^. 



2. Random partition of 11 and measurability 

In this section, we will consider random sets of points within 11 which essentially 
determine the process Z. Separating these critical points of 11 from the other ones, we 
get a random partition of 11. We will show that this partition is measurable, which 
allows for further investigation of this partition. 

For some fixed {t, z) e R'' x (0, oo), define the set 

Kuz = ({x,yJ)eS: f{t-x)>0, y= . 1, 
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which we call the set of points generating (t, z) due to the fact that 



z{t)^z ^ |nnX(,,| > 1 A nni^t,, = 0. 

Here, K = [jix.y.f)eK{^} ^ iv^ ^ if} ^or a set K C S. 

In a next step, we consider n fixed points {ti, zi), . . . , z„) S R'' x (0, oo). For 
any vector t = (ti, . . . ,tn) S K'* " and any mapping g with domain dom{g) C M.'^ we 
will write g{t) instead of {g{ti), . . . ,g{tn)), for short. Similarly, t > is understood as 
ti > 0, i = 1, . . . ,n. 

Now, we define the set of points generating (t,z) as 



K, 



\ {x, yJ)^S: max f{ti - x) > 0, y = mh 

I i—l,....n 



■ ,n f{ti - x) ) 

= e S -. yf{t - x) < z, yf{tj -x) = Zj for some j G {1, . . . 

This implies 

Kt,z. |(a;, y, /) e 5" : yf{tj - x) > z^ for some j € {1, . . . , n}| 

and 

^4:t,z n Xt.,^, = {(a;, yJ)eS: yf(t, - x) ^ z„ j//(t - x) < z} . 
Therefore, we have 



z{t) < z n n i^t.z = 

and 



Z{t)=z \UnKt,^,^^nKt,^\>l,i^l,...,nAUr\Kt.z = ?). (2.1) 

Now we define a random partition of 11 by 



Hi 
and lis 



= n n Xt,z(t) 
= nni^t,z(t), 
= n\(niun2). 



Relation (|2.ip implies that Hi = and n Ki.. z{ti)\ ^ 1 a.s. for z G {1, . . . , n}. 
Note that the partition of 11 into Hi, 112 and II3 corresponds to the classification of 
and in p . 

Before proceeding any further, we need to proof that Hi, 112 and II3 are well-defined. 
We will do this by showing the measurability of a further partition of 112 , namely the 
restriction of 11 to intersection sets. For any A e 2^^' - '"'^ \ 0, these are defined as 

Ia{z) = /^t.z n fl ift,,,. \ U Kt^,,^ 
\ieA jeA" 

= {{x,yj) G S : yf{t^ - a;) = z,, j G A, yfitj ~ x) < Zj, j ^ A\ 
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By construction Kt.z is a disjoint union of Ia{z), A <E 2^^' - '"^ \ 0. 

To proof the measurability of these restrictions of 11 to the intersection sets, let C 
be the a-algebra on generated by the cyhnder sets 

where si, . . . , s,„ e R'^, B e B"\ m e N. 
Proposition 2.1. Let ti, . . . , t„ E R'^ be fixed. 

1. The mapping 

is (B'^ X {Bn (0, cx))) X 2'^ ,C) -measurable. 

2. Let A e 2{i-- '"> \ (!) and B C S a bounded Borel set. Then, |n n LA{Z{t)) n B\ 
is a random variable. 

3. Hi, 112 and 113 are point processes. 

Proof. 1. It suffices to verify that 5'^^(Csi,...^s„(X™ ^^(ai, 6^))) is measurable for 
any Sj e R'^, Oj < 6^ S R, j = 1, . . . , m, m G N. We have 

bi \ 

X {/} 



U U (.-r^™7?r 



u u 

gi<g2 



/(t, -t)\=l.....m/(t,-t) 

(91,92) C ( : 



X (gi,'72) X {/}. 



As each / e G is measurable, sets of the type {t e R'' : f{ti — <) G -B} are 
measurable for any B G 6. Therefore, *"^(C,.i,....s„(X™ ^(ai, 6,))) G S'^ x (6 n 
(0,00)) X 2'='. 

2. We consider 

{w : |nn/A(z(t)) nB| = A:} 

= u n f u (-^"^ z{t)ExuU-i 



noGNrn— no yG 



,2/i H 

TO TO 



nn 



TO TO 



i G A, 



f{ti)<y,~^. 3iA\]c^B 



k 



By the first part of this proposition, is a measurable mapping and we get that 
{w : |n n LA{Z{t)) n = fc} is measurable. 

3. For any bounded Borel set B C S the second part of this proposition yields 

that iBi n s| = 0, |n2 n s| = EAe2^^ '.}\0 |n n /^(^(t)) n b\ and jng n b\ = 

|nn-B| — r\B\ are measurable. Thus, Hi, 112 and II3 are point processes ([1], 
Cor. 6.1.IV). 
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3. Blurred sets, scenarios and limit considerations 



This section mainly deals with the analysis of the distribution of the set of critical 
points, 112. As this set has the intensity measure zero conditional on an event Z{t) = z 
of probability zero, this distribution cannot be calculated straightforward. We need 
to borrow arguments from martingale theory, taking limits of probabilities conditional 
on the observations being in small intervals containing z. By this condition, the set 
of critical points gets blurred. We distinguish between different scenarios denoting 
which points influence the different observations. Using general bounds for the rate 
of convergence of the intensity of these sets, we can prove that each observation is 
generated by exactly one point of 11 (Proposition 13. ip . which restricts the number 
of scenarios that occur with positive probability. According to the blurred sets the 
scenarios get blurred, as well. The blurred scenarios are not exactly the same as the 
scenarios conditional on blurred observations, but much more tractable. However, it 
can be shown that both events asymptotically yield the same conditional probability 
(Theorem 13. ip . Based on these considerations, the independence of Hi, 112 and lis 
conditional on Z{t) is shown (Theorem 13. 2p . This allows us to simulate 112 and lis 
independently. Furthermore, Ha turns out to be easily simulated (Theorem [ 



Let 



-Tvn. — O 



Z{u) e — ,— 



Then, {Fm\m&i is a filtration and Foo 
jm{z) e N be such that z G Am{z) with 



,n, k eNo 



^ ^\meN^r,^ = <y{Z{t)). For z > 0, let 



' \ 2m ' 2™ 



Then, we have Amiz) ™ — f {z} monotonically and {cv £ fl : Z{t) G Am{z)} G J-„ 
Furthermore, let 



and Kt. 



-(m) 



U Kt^, = {{x, yJ)eS: yf{t - x) & A„(z)}, 

>„(z) + 1 



U ^M-{(x,y,/)G5: y/(t-x)< 



yf{U - x) e Am{zi) for some i G {1, . . . 
{x,yJ)eS:yfiU-x)> ^"'^':l^^ 



for some i G {1, . . . , n} > = X j 



These definitions imply that 



-(m) 



and kI;^^ U Xt,z"'" - 
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We call these sets the blurred sets belonging to Z{t) conditional on Z{t) E A„i{z) 
This notation is due to the fact that we have 

>^,(z) + l ^ 

Furthermore, as 

we get that 

Thus, we obtain 

z(t)eA„(z) ^ |nn<™i^nifi™^|>i, z = i,...,n a nni?;:;^™^ - 0. (3.1) 

In particular, for fixed z G (0, cxd)", the point process n\(ift';^^Uift.z ) is independent 
of the event Z(t) G A,„(z). 

Based on these blurred sets, we define the blurred intersection sets 

(z) = if n n \ U A e 2{i.--"> \ 0. 

We note that i^t™' can be written as a disjoint union 

Lemma 3.1. For any A E 2{i- ""> \ and z > we /laue A (^/^"'(z)j e 0(2-™), 
where A(-) is given by (|1.2p . 

Proof. It suffices to show that 

A(ifl"i^) eO(2-™), » = l,...,n. (3.2) 
By a straightforward computation we get 



1 



^f( F{ti - x) dx 
\Jr'' 



1 

2'" 



, 2™ 2™ / 2™ 2" 



Using the fact that lini,„_>.oo '''^m''* = Zi, the assertion of the Lemma follows. 

In Section m a more precise notion about the speed of convergence of im{zi) 
Zi will be useful. 
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Lemma 3.2. For any e > 0, with probability one we have 

liminf2'"(i+-) min ( il^li^Mll ^ z (ti) ] ^ oo 

m— >-oo i— l,...,n \ 2^ 

and linisup2'"(i+^) max f ^I^ii^M _ ^(t, 

m— >oo i— ±,...,u y ^ 

Proof. We present the proof of the first assertion. The second one can be shown 
analogously. For the first assertion it suffices to show that 

liminf 2™(i+^) f ^^^i^^M±l _ zit,)) = oo, ^ = l,...,n. 
Let a G (O, I) and m G N large enough such that 

l^-^p^^^p (-^2^) ^ i^^^ (-^) ^ ' ^'-'^ 

Then, we have 

= ^ P ((/s - a)2-" < < /c2-") 

1 \ / / I a/fc 



< 



fc=i 

oo 

1 \ / 1 a 



k=l 



where we used the fact that 1 — exp(— a;) < x for all a; > 0. Employing p.3p . we get 
for a = C2~™^ with C > and m large enough that 

P (^llA^Ml±l _ z(tO < C • 2-™(i+^)^ < 4C • (2^)-™. 

Therefore, the probabilities above are summable with respect to m and the Borel- 
Cantelli lemma yields 

P ( liminf 2™(i+-) ( + 1 _ z{t,)) <C)^0 

\ m— >oo \ 2™ / / 

for any C > 0. This completes the proof. 



Now, we relate the points in K^"^' to the corresponding "blurred" intersection sets 
by introducing disjoint "blurred" scenarios 

= |nn/i'"'(z)| =71^, Ae2{i'-'"}\0, |nni?^^"^| = 0} 
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with {uA, A £ 2{i' -"> \ 0} e TVi where 



A^i - I {riA, A e 2{i-"> \ 0}, ^ > 1, * = 1, 

I A-.i&A 



Thus, we have 



{Z(t) e A„(z)} = y £;f:]j(z) 



and the union is obviously disjoint. In the same way, dropping the (m) in the definition, 
we specify scenarios i?{„^}(z). 

Now, we show that 

P(3*e{l,...,n}:|nni^,,,z(t,)l>2)-0. 

To this end, we first verify that, with probability one, 112 has no points which can be 
removed without any effect on Z{t). We consider scenarios i?{„^} (2'(t)) with {ua} € 
N2. Here, N2 is the ensemble of all sets {riA, A g 2^^'- -^"'^ \ 0} such that there exists 
some A* G 2{i' - '"> \ with 

^ inA'-l, A^A*, 
^ ytiA, else. 

and \n*A\ G ^i- Thus, {n^} G if and only if i?{„^}(Z(t)) remains an allowable 
scenario after removing one point. 

Lemma 3.3. With probability one, we have 



p I u i?{„,}(z(t)) 



lim P 11 E^r\(z) 

.{nA}eN2 



Z{t) = z 

z(t) e ^„,(z) I = 0. 



In particular, hm,„^oo P (|n D if^^J^ \ Hal > Z{t) E A,„(z)^ = 0. 

Proof. By Levy's "Upward" Theorem (see [TB], Thm. 50.3), we a.s. have 

^ P (Z(t)) I Z(t)=z)=Jim^ ^ P(i?{„,}(Z(t)) |Z(t)eA„,(z)). 

Noting that 

u u 

it suffices to show the equation 



hm P(i?:i}(z)|^(t)eA„,(z))= hm 



{nA}eN2 



0. (3.4) 
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for every z > 0. Corresponding to each i?|™^j(z) with {ua} G N2 there exists some 
£;|™ij(z) with {n*jj e A^i as in the definition of N2. For A* e 2{i'- \ let N2,A' 
be the set of all {ua} G ^2 with n\, = tla* — 1- Then, for fixed to e N, the left-hand 
side of p.4p is less than or equal to 

(m) 



E 



p^:i}(-)) 



A-62{i.->"}\0 E{ n A } e Wl 



E 



A.e2U.-.-.">\0 ^ ,-A(/r) A(/.Jz))-. .exp(-A(<"))) 

{«A}eAri \_Ae2{i "}\0 / 



< 



E 



{ri*i}eAri Ae2{i.--."}\0 



A-e2{i.-.'.}\0 E n ""nV 

{KAleATi Ae2{i '•}\0 

< A(/i?)(z)), 

A*G2{1.' ■."}\0 

where we imbedded N2^A' hi^to A^i by identifying {n^} with the corresponding {n^}. 
By Lemma [5TT1 we get equality p.4p . i.e. the first assertion of the lemma. Furthermore, 
Inni^t™' \n2| > O and Z{t) G A„i(z) imply that IlnK^^^ contains points which can 
be removed without affecting Z{t) e A,n{z). This is, 

lim p (|nn ii-t^™^ \ Hal > I Z(t) e A™(z)) 

^jl- E P(i?E}(z)|^(t)GA„(z))=0. 

{nA}<£N2 

Thus, the second assertion of this lemma is verified. 

The following proposition is also stated in a more general setting in T, Prop. 2.2. 
Proposition 3.1. For any fixed ii, . . . , t„ e R'' we have 

P(|n n Kt^,z{u)\ > 2 for some i E {1, . . . ,n}) = 0. 

Proof It suffices to show P(|n n Ku,z{u) \ > 2) = for ah i e {1, . . . , n}. Now, let 
i e {1, . . . , n} be fixed. Then, conditioning on Z{ti) only, we get 



P{\IlnKt^,z(u)\>2\Z{U)^z)^F\ y E[,,^}{Z{U)) 



Z{U) = z I = 

for almost every z > by Lemma 13.31 This yields the desired result. 
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Figure 1: The event Am in the proof of Thm.[3TT] Black crosses: data {ti, Z{ti)), black line: 
Kt,z(t), grey area: K^'^zit)^ black dots: H n Tf^'^'j^j. Left: |n n A'^^'^'t, \ > 0. Right: 

|nn/{i}(^(t))| = i>o = |nnj|rAz(t))| 




0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 



Proposition 13. ll ensures that almost surely one of the scenarios i?{„^}(z) with 
Sa- i<£A "-^ = 1 foi' alH € {1, . . . , n} occurs. 

Theorem 3.1. With probability 1 we have 

1. P(i?{„,}(Z(t)) I Z{t) = z) = Jim^P(i?|7]j(z) I Z{t) e A„,{z)) 
for any {ha} G Ni, 

2. P (|n2 n rj, J = 1, . . . , fc, iHa n B,\ =r,, j = k + l,...j\ Z(t) = z) 

= lim p (|n n K^":^ n r„ J = 1, . . . , fc, 

|n n B, \ u = r,, j = fc + 1, . . . , / I z(t) e A„(z)) 

/or any C S, r^ G N, j G {1, . . . , /}. 
Proof. 1. Again, Levy's "Upward" Theorem ([16], Thm. 50.3) yields 

P{E^^,}{Z{t)) I Z(t) - z) - lim P(£;{„,}(Z(t)) I Z(t) e A„(z)) 

with probability 1. It remains to verify 

lim PiE'^^Uzit)) I Z{t) e A„(z)) - lim P(i?{„,}(Z(t)) | Z(t) G A™(z)) 

(3.5) 

where P(£;|;:]j (Z(t)) | Z{t) € A„(z)) equals P(i?|r]}(^) I Z{t) G A„(z)) by 
definition. 

To this end, we consider the symmetric difference A^, of the events £'|™j|j (Z(t)) 

and i?{„^} (Z(t)). Note that any element of A,„ satisfies |n H ^'^["^(t)) \ ^2! > 

(cf. Figure m left) or |nn/A(^(t))| > |nn /^"^(Z(t))| for some A e 2{i-'"> \0 
(cf. Figure [U right). The second kind of event happens if there is a point of 11 
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in /^(Z(t)) n (Uji^A ^t"z(t ))■ t'^i^ vanishes for any Z{t) > as m — > cx), 
we get that 

{a;ef7: |n n /^(^(t))! > |n n /^"\z(t))|} \ 0, m^cx), 
for any A € 2^^^- -^"^ \ 0, and therefore 

hm P(|nn/A(z(t))| > |nn/i'")(z(t))|) = 0. 

m— )-oo 

This yields 

/•oo 

/ hm p(|nn/^(z(t))| > |nn/i"^(z)| | z(t) e A™(z)) P(z(t) e dz) 
= hm p(|nn/A(z(t))| > |nn/i"^(z(t))|) =0 

using dominated convergence and the fact that 

/•oo 

p(|nn/^(z(t))| > |nn/i")(z)| | z(t) e An(z)) P(z(t) e dz) 
= ^ p(|nn/^(z(t))| >|nn/i")(z(t))| I z(t)e A„(z))p(z(t)e A„(z)) 

= p(|nn/^(z(t))|>|nn/i")(z(t))|). 

Therefore, we have 

hm p(|nn/A(z(t))| > |nn/^"^(z)| I z(t) e A„(z)) = o (3.6) 

for any A g 2^^' "'"^ \ and almost all z > 0. All in all, we end up with 



lim 

>-00 



(4:]j(z(t)) I Z{t) e A,„(z)) - p (i?{„,}(z(t)) I Z{t) e A,„(z)) 



< hm V P (\iiniA{z{t))\ > |nn/i")(z)| | z{t) e A„,(z) 

+ lim P (|n n \ Hsl > I Z(t) e A„(z)) = O a.s. 

by (|3.6p and by the second part of Lemma [5?^ Thus, we get p.Sp which completes 
the proof. 

2. Let Bi, . . . , Bk,Bk+i, ...,BieB'^x{Bn{0, ooj) x 2"^. Then, each of the events 

{|n2 n BjI ^ |n n K^''1\^^ n S^l for any j = l, . . . , k} 
and {|n3 n ^ |n n \ (Xg}^) U K^)''""^)\ for any J fc + 1, . . . , 
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implies that |nnX^'^|^.^ \n2| > 0. Therefore, for any ri, . . . , r^, Vk+i, . . . , r/ G N, 
we have 

hm p(\U2nB,\=r,, j = l,...,k, 

iHa n Sjl = r,, j = fc + 1, . . . , ? I Z{t) e An(z) 

-p(|nnx;,'"j)ni3,Hr„ 

|n n Bj \ (i^g) u = r,, J = fc + 1, . . . , z I z(t) e A„(z) 

< hm p (|n n \ n2| > z(t) e A„,{z)) = o. 

by the second part of Lemma 13.31 This verifies the assertion. 

These results enable us to show the independence of Hi, 112 and II3 conditional on 
Z(t) and to calculate the conditional distribution of II3. Similar results can be found 
in [7], Thm. 3.1, obtained by a different approach. 

Theorem 3.2. 1. With probability 1 the point processes T\i, 1^2 andTi^ conditional 
on Z{t) are stochastically independent. 



2. The process \ Z{t) — z has the same distribution as 11 \ {K^ .^, U Kt^z) with 
probability 1. 

Proof. 1. Since Hi = a.s., we only have to show the independence of 112 and 
conditional on Z{t). Let Bi, . . . , Bk, Bk+i, . . . , Bi e B"^ x ((0,cx)) D B) x g 
and ri , . . . , r/ G Nq . The second part of Theorem 13.11 yields 



P(|n2ns,| = rj, j = i,...,fc, in-^n Bj\ = rj, j = k + i,...,i\ z(t) = z) 
= jirn^ P(|n n ifg) n r„ j = i, . . . , fc, 

|n n B, \ (xg) u x^^"^)! ^r,, j = k + i,...,i \ z{t) e A™(z)). 

By dSH) the process U \ (K^'^^ U X^^™') is independent of the event Z{t) e 
A,„(z). Hence, 

p(|n2nB,| =r,, j = i,...,fc, ingns.Hr,, j = k + i,...j \ z(t) = z) 
= lim P(|n n n = rj, J = 1, . . . , fc I z(t) e A„(z)) 

• lim P(|n n B, \ u = r,, j - fc + 1, . . . , z I Z(t) e An(z)) 

= P(|n2 n Bj\ = r,, J - 1, . . . ,fc I z(t) - z) 

• P(|n3 n B,\ =r,, j = k + i,...,i \ z{t) = z). 



where we use the same arguments as before. 
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2. For any sets Bi, . . . , Bi £ B"^ x {B Ci {0, oo)) x G and ri, . . . , n £ No, the second 
part of Theorem 13.11 impHes 

P(|n3nBj|<r„ j = i,...,i\ z{t)=z) 

= hm P(|n n B, \ {K^^l^ U < r,, J - 1, . . . , ? I Z(t) e A„,(z)) 

771— >-00 : I 

= hm P(|n n B, \ {k'^^ U < r,, J = 1, . . . , 



= P(|nnB,\(ift,zUift,z)| <rj, j- = i,...,o. 

Here, we use that the process 11 \ {K^^^ U i^t.z^™'') is independent of the event 
Z(t) E Am{z) and the fact that we have a monotone hmit. This completes the 
proof. 

By the second part of Theorem [32] the process lis | Z{t) — z can be easily simulated 
by unconditionally simulating 11 and restricting it to M'' x (0, oo) x G \ (i^t,z U A't.z). 

At the end of this section we note that there exists a more general version of Theorem 
IS.ll which we will need for simulation. Let Bi, . . . , Bk, G S'^x ((0, oo)nS) xC/ be pairwise 
disjoint with — We introduce generalized "blurred" scenarios 



Fy ' (7) - Fy ' (7) 
= {|nn/f'(z)nB,|=n^\ Ae2{i-"}\0^ j = i,...,fc, inniT^^™^! = 0} 



with n^J^' e No such that Y. -=i Ea: teA "a > 1 for i G {1, . . . , n}. 

Analogously, generalized scenarios £'|^(j)j(z) are defined. Then, in exactly the same 
way as Theorem 13.11 the following theorem can be shown. 



Theorem 3.3. With probability 1 we have 

P(i? „u, (Z(t)) I Z{t) = z) = hm P(i?;";)„, (z) I Z{t) e A^{z)) 

for any scenario £'^^o)j(z) with Ej=i J2a: isA "a^ ^ 1 /o'^ all i £ {1, ... , n}. 

The remainder of the paper will address the problem of simulating 112 | 2^(t) = z. We 
propose a procedure consisting of two steps. First, we draw a scenario £'{„^}(Z(t)) 
conditional on Z{t) = z. Then, the points of 112 corresponding to this scenario are 
simulated. 



4. Calculations in the case of a finite number of shape functions 

As shown in Section [3l all we need for calculating P(i?{„^} (Z(t)) | Z{t) — z) is 
the exact asymptotic behaviour of A{I^\z)). In particular, we have to analyze the 
behaviour of the intersection of two curves Kt-^zi+Si H Kt-^zj+Sj for |(5i|, \6j\ small. 
Explicit calculations turn out to be quite involved. Therefore, we restrict ourselves to 
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Figure 2: Blurred intersection sets for \A\ = 1, |j4| — 2 and \A\ = 3. Black crosses: data 
(ti,Z{ti)), dashed black lines: -K^t •,z(ti) , black line: Kt.z{t), grey area: K^"^^^y black area: 
7^'"^(Z(t)) with \A\ = 1 (left), \A\ = 2 (middle) and \A\ = 3 (right). 




0.0 0.5 1.0 1.5 2.0 00 0.5 LO 1.5 2.0 OO 05 LO 1.5 2.0 



the case d = 1. Furthermore, the intersection depends on the derivative of the shape 
function. We calculate the asymptotics of A{l'^\z)) for \A\ — 1 (Proposition 
1^1 = 2 (Proposition WA} and |^| > 3 (Proposition IT^ . see Figure [5] In the latter 
case, the rate of convergence of A(/^™'') cannot be determined exactly. Nevertheless, 



the conditional probability of any scenario can be calculated (Theorem 14. II) . 

First we assume that G is a finite space of functions / : M — > (0, oo) such that the 
intersections 

Me,t„ ={teR: f{t) = cfito + t)} (4.1) 

are finite for all c> 0,to eR, f e G. This implies that each set /a(z), A e 2{i--"} \0, 
\A\ > 2, z > 0, is finite. W.l.o.g. we assume that Pi?({/}) > for all / G G. 

Proposition 4.1. Let ti, t2 G R, zi, 22 > such that 

I{i,2}{2) = {{to,yoJ)}- 

Furthermore, let f be continuously dijjerentiable in a neighbourhood ofti — to andt2 — to 
with 

zifit2-to)^Z2fit,-to). (4.2) 

Then, we have 

9 — 2m 

A(/|;1}(z)) - J- PF{{f}) + 0(2-2"). 

^ ■ ^ yo\zif'(t2 - to) - Z2.f'iti - to)\ 

Proof. We note that {tg, yo, f) satisfies the equation 

fih - to) .f{t2 - to) _i 
= = Vo ■ 

Zl Z2 

Let 

iJ : (-z,oo) xM^R, {S,t)^ /(^i " ^) _ 

Zl +O1 Z2+ 02 
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Then, H{0,to) = and 

— (0,to) = --^ + + 

at zi Z2 

due to (I4.2p . The impHcit function theorem yields the existence of a neighbourhood 
y of and a continuously differentiable function /i : 1/ — > M such that H{5, h{5)) — 0. 
Using the notation {ts,ys, f) = I{i,2} (z + S) we get h{5) — tg and the equality 

f{h - ts) ^ f{t2 - ts) ^ _i 

Zi +Si Z2+ $2 

As h is continuously differentiable, we obtain tQ — ts G 0{\\d\\), and a Taylor expansion 
of / yields 

fiU - ts) = fit, - to) - fit, - to) ■ its - to) + o{\\6\\). (4.3) 
Let git) — fit — to). Thus, using (|4.3I) . ts is given implicitly by 

giti)-g'iti)-its-to) ^ git2)-g'it2)-its-to) ^ 

Zl + Si Z2+ S2 

which implies the explicit representation 

+ + '^1.9(^2) - S2giti) 

ts = to + y— —— + oi\\5\\). (4.4) 

zig it2) - Z2g'iti) 

Plugging in (lO)) into yields 

-1 fih-ts) giti) g'iti) 6igit2) - S2giti) 

ys " = — rii T7r\ T7T^ + oi\\6\\). 

zi + di zi+di Zl zig'it2) - Z2g'iti) 

As / and 6 1-^ ts — /i((5) are C^-functions, all the terms o(||(5||) are continuously 
differentiable for small \ \S\\. Therefore, the mapping 

$:F^Mx (0,00), S^its,ys^) 

is continuously differentiable near the origin. Calculating the partial derivatives ex- 
plicitly we obtain 

detiOm) = - 2 I f,t^ 77T^ + (4-5) 

zt ■ (219 (^2) - Z2g'iti)) 

As det(Z?$(0)) 7^ 0, the inverse function theorem allows to regard $ as a diffeo- 
morphism restricted to an appropriate neighbourhood of 0. Thus, considering the 
transformed Poisson point process 11 = u)en '^(s,«-i) on M x (0, 00) whose intensity 
measure is the Lebesgue measure, we get 

A({/{i,2}(z + (5), Zi + S,eA,niz,), i = l,2}) 

Pf({/}) dit,y) 
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\det{D^S))\-FFi{f})dd 



(Ar^(zi)-^zi)x{A^(z2)~Z2) l^lS'l^s) " 22.9'(^l) 



PF{{f})dS. 



We note that the term o(l) is contmuous w.r.t. S and therefore the integrands can be 
locally bounded by 



\z,g'{t2)~Z2g'ih)\ 



from below, and by 

Wo 



\z,g'{t2)-Z2g'{h)\ 



Pf({/}) 



Pf({/}) 



from above for all (61,62) G {Am{zi) — zi) x (Am (22) — Z2) with m large enough and 
an appropriate sequence {Sm)m<£N with Sm \i 0. This implies that the integral has the 
form 

\zig'{t2) - Z25'(ii)| 
which is the desired result. 



Remark 4.1. 1. Using j:^^^^^^^ — f(^t2-to) ^ ^0 we get that the equality 



Zlf{t2-to) = Z2f{h-to) 



holds if and only if 



d z\ 



dt f{ti - t) 



d Z2 



dt f{t2 - t) 



i.e. if and only if the two sets of admissible points, i^t^.^^ and Kt2,z2 1 ^-re tangents 
to each other in {to,yo,f) which is an event of probability zero by Assumption 
(|4T|) . Therefore, (|t!2l) is satisfied a.s. 



2- If -^{1.2} (2;) consists of a finite number of points, 

/{i,2}(^i, ^2) = yo'^ ■ • • : (4", yo'\ A)}, 

we get 

' P({/.}) 

,=1 %'-i^i/i(i2-4'V^2/i(ti-4'^)i 



+ 0(2-2™). 



Proposition 4.2. Let . . . , G M, zi, . . . , > 0, / > 3 such that 

-^{i,...,;}(z) = {(io,yo,/)}- 
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Let / : M — > (0, oo) be continuously dijferentiable in a neighbourhood ofti — tg, . . . , U—to 
with (j4.2p . Then, we have 

For any C > 0, e > 0, there exists mc.e G N such that 

A(/|™)^ _;j(z)) >C2-2'»(i+^) 

for all m > mc.e- 

Proof. The first assertion follows immediately from Proposition |4?1] by the fact that 

I 

tA™-) r- tA^^) Z. 



In order to verify the second assertion, we recall results from the proof of Proposition 
14.11 we showed the existence of a C^-function (5i , ^ t5iS2 defined in an appropriate 
neighbourhood V of (0, 0) such that 

fjti ~ tsiS2) _ I{h - ts^Si) 

Zi + 5i Z2+ S2 

Now, we consider the C^-functions 

zi + Oi Zi + di 

As iJ,(0, 0, 0) = and |f-(0, 0, 0) = £lh^ ^ 0, we get the existence of a continuously 
differentiable function hi defined on a neighbourhood of (0, 0) such that 

fjtl - ts,S2) ^ fjtj - ts-^Si) ,^ gx 

zi+5i Zi + hi{Si,62)' 

Using Taylor expansions of g{-) = /(• — to) of first order, employing Equation (|4.4p . 
and solving Equation (|4.6p yields 

h fx X \ 9(M)c . z,g'{ti) - zig'{ti) gfta)^! - g{ti)52 . , 

hi[di,d2) = -7rT''i + 771 7771 TTn + o(\0i\) + o(\62\). 4.7) 

5(^1) 9[ti) zig'{t2) - Z2g'{ti) 

So, there are constants ci.i, C2.i such that hi{6i, 82) = ci.iSi + C2.i52 + od^^il) + o(|<^2|)- 
Let Am = Am{zi) — Zi for z e {1, . . . ,ri}. We are interested in those pairs {5i,52) G 
X Aln^ with h,{Si,S2) e An-*. By Lemma [321 for any C" > 0, e > 0, we have 
that (-C"2-"(i+'^),C"2-™(i+^)) e A^ri}, i = 1, . . . , n, for m large enough. Therefore, 
hi{Si,S2) e Am is guaranteed for \Si\ < ^ \\c^ ^ind \62\ < ^ \\c2 -l^^ 
sufficiently large. 

By the same argumentation for all i e {3, ...,?} we get that the existence of all 
hi{Si, 62) is ensured for 

\^^\<T^ \S2\<^ 1 1 (4-8) 

3maxi^3 7 \ciA Smax^^g , \c2a\ 
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for m large enough. Furthermore, to ensure 5i £ Am\ 62 £ Am\ we have to add the 
conditions \S1l\S2\ < C"2~"(i+'^). With Cj = max{l, 3maxi=3,...,/ |cj- ,|} for j = 1,2, 
this yields 

A({/{i,....,}(z + <5), <5, ev4W, z^i,...,/}) 

yiCiC2\zif'{h-to)-z2f'{h-to)\ ^ ' 

where we use the same argumentation as in the proof of Proposition l4.1l This completes 
the proof. 



Remark 4.2. We note that Proposition l4.2l does not provide exact asymptotics which 
we need to determine P{{to,yo,f ) e 11 | |nn /{i,...,/}(z)| = l,Z{t) = z) for (to,2/o,/) e 
s-iid |/{i_..._/}(z)| > 1. However, we can get exact results by conditioning on 
Z{ti) being in intervals of different size for each i G {1, . . . ,1} instead of Z{ti) e Am{zi) 
for all i — We will choose these intervals such that some restrictions on 

the intersection sets vanish asymptotically and we can resort to the results on the 
intersection of two curves. 

The calculations in the proof of Proposition 14.21 vield 

M6,,52)\ < (|C1,,| + 0(1)) -1(511 + {\C2.,\ + 0(1)) • |^2| < 2-"(|ci,,| + |C2,.| + o(l)) 

for {61,62) S Am^ X Am'. Thus, for any e > 0, using the same arguments as in 
the proof of Lemma [3.21 we can replace m by [m(l — e)J in Equation (|4.8p and get 

that h,i6i,62) G ^olds for \S,\ < ""''"ticjr'^" ^ 2-'™2-" and \62\ < 

^ ^ '^3|c2 1'^' ^ ' ~ ^"^ rn2-in fgj- large enough. Therefore, 

/ii((5i,(52) C i = 3,...,l 

for all 61 G A^^' C (-2-'", 2-™], 62 G aI^^ C (-2-'",2-™] if m is sufficiently large. 
This implies 

+ -JigAW, 62eA<^J, '5,eA«(^_^)j, i = 3,...,/} 
= {/{i,...,}(z + (5), (5i,(52) G AW X A(^), /i,(^i,<52)G A«(i_,)j, * = 3,...,/} 

and, therefore 

A ({/JS.,,}(z + 5), ,5i G A,W, ^2 G A(2), 6. e A« , z = 3, . . . , z}) 

2--PH{/}) ,^(2_2m) 



y§ki/'fe-<o)-^2/'(ii-io)l 
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By conditioning on Z{t) e Ajn{zi) x Arn{z2) x Xi=3,...jylL„(i_e)j (z^), for /{i,...,/}(z) = 

{(4^^2/o^^/l)>•■•)(4''^yo''^/fc)}' ^ > 3, we can apply Levy's "Upward" Theorem 
f|16). Thm. 50.3) and end up with 

P((io,2/o,/) G n I |nn/{i,...,i}(z(t))| = i,z(t) = z) 

|nn/{i,...j}(z(t))| = 1, z{t) e A„(zi) x A„(z2) x x,=3,...,iA Lm(l-e)J (^i) ) 



Pf({/}) Pi.({/,}) 



(4.9) 



Note that Levy's "Upward" Theorem imphes that, with probability one, the right- 
hand side of (|4.9p does not depend on the choice of labehing. 

Thus, despite the lack of exact convergence rates of A(j|™'' ;j(z), the distribution 
of nn J{i^...,i}(2'(t)) I |nn /{i,...,;}(Z(t))| = lean be determined exactly. 

Proposition 4.3. Let / e G, t e M", z > such that f is continuously dijferentiable 
in a neighbourhood ofti—tQ for alii £ {1^ ... ,n} and all {t^, ya, f) G Kt.z(\j?=i{Kti .zi^^ 

Kti^zi), i-S- cili (^Oj2/0:/) € Kt-^zi involved in an intersection. 

We denote the projection of the set /{,j}(z)n(Rx(0,oo)x{/}) onto its first component 
in R by 

D\f^ = {< e M : (t, y, /) e /{,} (z) for some y > 0}. 

Then, we have 



A 



(j;j(z) n (R X (0,00) X {/})) = 2-".P^({/}). / l^^^ At+o{2-^). (4.10) 



Proof. First, we note that by renumbering it suffices to show the result for i = 1. 

The idea of this proof is to assess the set Dl by the sets D\ J^^^ from below and 
^i™nax from above. Here, d["^-^ consists of all first components of l'^'^ {z) which 
are not part of any intersections l'^\z), A D {1} and i?i™]ax °f ^^^^ 

components of 1Jad{i} Analogously, A ^/|"j^(z) n (R x (0, 00) x {f})j can 

be bounded from below and above by replacing d[^^ in (|4.10p by the sets d["^^-^^ 
and i'l'^ax 1 respectively. We will show that the difference, which consists of blurred 
intersections I^\z), A D {1}, vanishes asymptotically. 

Let Am = Am{zi) — Zi. Then, for any S G X"^]^(— z^, 00) we define 

n-]^{te.: (,^,/).W^ + ^^)} 

tem: ^±^< min . (4.11) 

/(tl-t) ^=2^...^nf{t,-t)\ ^ ' 
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Thus, defining d'(%^ = n,ex?^,A« ^ul^ ^ U,ex?^,A« ^M' 

On the other hand, we have 

[{t,vJ) e R X (0,oo) X {/} : t e d[%^^, yfih - t) e A^izi)} 

'{1} " 



C l[T.\z) c{{t,yJ)eRx (0,oo) X {/} : te D^/X^^, yf{h ^ t) e A„,(zi)} 



(4.12) 

Now, let t £ D^^l^^ \ . Then, by definition of d[X.^ and D[fl^^^ , there exist 

,5(1), (5(2) e xf^^A^^' such that t e but < (f. D[^li2)- That is, by Equation (|iTl|) . 

— ■ < mm — —— and — ^ > mm — -—. 

fih - t) t=2,...,n f{U - t) fih -t) - ^=2,...,n f {t, ~ t) 

By continuity arguments, a (5 G x"^j^Am exists such that 

zi+ 5i . Zi + 5i 

mm 



f{h-t) ^=2,...,n J{U-ty 

i.e. t e T^"^ = e M : {t,yj) e Ua: {^gaIa^H^) for some y > 0}. Thus, 

<Ja.\<LcT('"). (4.13) 

By definition, J!*-™) denotes the set of first components involved in any blurred inter- 
section and we have T^('") \ Ti = {t e M : (i, y, f) e IJ^. {ijcA -^'4(^) fo'^ some y > 0} 
as m — >■ cxD and Ti is finite by Assumption (j4.ip . Therefore, dominated convergence 
yields 

f{ti-t)dt\Q, m~^oo. (4.14) 



Thus, by Equations (|4.12p and (|4.13p we get 
a(/{7j)(z)A {(t,y,/)eMx(0,c^)x{/}: t G 2//(ti - t) G A„(zi)}) 

<A({(i,2/,/)eMx(0,(x))x{/}: te^^Lx, ^/(ti - e ^m(^i)} \ 

{(t, y, /) e M X (0, c^) X {/} : i e i^^^^n , vfih - e An(^i)} ) 
< A ({(i,t/, /) e M X (0,oo) X {/} : te t['^\ yfih - t) e An(zi)}) 

The last equality follows from Equation (|4.14p . 
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Hence, we have 

= 2-".. PHm)-(/^„,^d, + „(!)) 

which completes the proof. 

Now, we can use Theorem 13. II in order to compute the conditional probabilities 
P(i?{„,}(Z(t)) I Z{t) = z) = Jim^P(i?|;:]j(z) I Z{t) e A„,{z)) 

^ lim ^ " (4.15) 

where No = {{ua : A e 2{i' -^"} \ 0} : J^^. = 1, z = 1, . . . , n}. 

As the sets /j["^(z), A £ 2^^' - '"^ \ 0, are pairwise disjoint and the sequences 
A{l'f\z)) tend to zero for TO oo by Lemma I3.1[ we get 

P(i;|:]j(z))^ exp(-A(i^)) n MlAiz)). (4.16) 

A: nA = l 

Considering (j4.15p . we can restrict ourselves to those scenarios with the slowest rate 
of convergence to zero. Propositions 14.11 14.21 and 14.31 yield that scenarios involving 
intersections of at least three sets are always of a dominating order. Therefore, the 
unknown terms from Proposition l4.2l are cancelled out. 

Example 4.1. Let F{x) = f{x) = --^exp(— ^) a.s. Furthermore, let Z{ti) = 
Z{t2) = a where t\ ~ —b and t2 = b for constants a, 6 > 0. Then, we get 

\2 ^ 



/{i}(a,a)= |(^<,a\/2^exp(^^-^^-^^ : i < o| 
/{2}(a,a)= I a\/27r exp 
I{i,2}{a,a)^ |^0,aV27rexp 



By the formulae from Propositions 14.11 and 14. 3[ we get 

A(/|^/(a,a)) = ^exp (^^^^j dt + o(2-™) ^ ^-^$(6) (4.18) 

2-2™(i 1 exp(-^))2+o(2-2") 2-2" 

A(/7,(a,a)) = ^ ^ ^= ^. (4.19) 

a^exp(-f)-a^exp(-f) 2a%V2^eM'-^ 
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Here, there are two scenarios satisfying \Ilr\Kt^^zi \ — \^<^Kt2,z2 1 — 1 and the Equations 
gT7D - (|iT^ yield 

p (inn = 1 I zi-b) . zit) ^ a) . ^ "^f ^^^^^^^.^^ ^ , 

P(|nn/{i}(a,a)| = |nn/{2}((a,a))| = 1 I 

^(-fe) = Z{b) =a)= . 

Note that the probabihty that the observations are generated by two different points 
of n increases as a gets smaUer. This is due to the fact that 11 gets more intense as 
the second component decreases. 

Using the formulae above, the limits of the conditional probabilities can always be 
calculated explicitly except for those cases where two scenarios exist, both involving 
different terms which cannot be determined exactly (cf. Proposition 14.21) . This may 
happen only if two sets Ai = {ii, . . . , ir} and A2 — {ji, . . . , js}, r, s > 3, Ai f] A2 ^ 0, 
exist such that 

Ja, {Zit)) ^% , Ja, {Z{t)) ^ and Ja.vja, {.Z{t)) = 0, (4.20) 

where 

JA{z{t)) = y iB{z{t)) = ift.z(t) n fl A e 2{i-^"> \ 0. 

BdA ieA 



In all other cases the terms as in Proposition l4.2l are cancelled out. Note that we work 
with sets of the type JA{Z{t)) in order to avoid case-by-case analysis for all the sets 
/s(Z(t)) with BdA. 

Lemma 4.1. Let G consist of Junctions which are continuously differentiable a.e. 
Then, for any fixed set {ti, . . . , i„} C M we have 

P{Z{t) satisfies = 0. 

Proof. We proof that condition (|4.20l) has probability for all fixed index sets 
Ai,A2 C By renumbering, we may assume Ai = {l,...,r} and A2 = 

{q, . . . , q + s — 1} with q < r. Assume that P{Z{t) satisfies (|4.20l) ') > 0. In a first step 
we only consider those realisations of Z{ti), . . . , Z{tr) with J^j(Z(ii), . . . , Z{tr)) 7^ 0. 
Then, by the calculations in Propositions 14. 1[ 14.21 and 14. 3[ we get that 

p(|n n /i7'(z(ti), . . . , z{tr))\ = 1) ^ o(2-2™(i+e)) 

for any e > and 

p(|n n 47^(z(ti), . . . , z{tr))\ = 1, |n n 4";H^(^i), ■ • • , z{tr))\ = i)e o{2-^n 

for any Bi,B2 C Ai, Bi D B2 ^ 0. This yields \U D JAA^iti), ■ • ■,Z{tr))\ 1 almost 
surely. 
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Similarly, for almost every Zq, . . . , Zq+s-i such that JA2{Z{tq), . . . , Z{tq+s-i)) 7^ 0, 
we have |n n Ja^ (^(ig), • ■ • , = 1- As 

{uj : Z{t) satisfies (g^Hl)} 
(1{uj: JA,(Z(ti),...,Z(t,)) ^0}n{w: JA,(Z(t5), . . . , Z(V.-i)) 7^ 0}, 

we have |n n J^^ (^(ti), • . . , ^(tr))| - 1 and \n D JA,{Z{tq), . ■ . , Z{tq+s-i))\ = 1 for 
Z{t) satisfying (I4.20p almost sm-ely. Therefore, we get P{\Iln K(^t.^z{ti)) \ > 2) > for 
every i £ Ai Ci A2 since J^i^uAa (^(t)) = 0- This is a contradiction to Proposition 13.11 
which completes the proof. 

From the considerations above and Lemma |4. II we immediately derive the following 
result. 

Theorem 4.1. Let G consist of functions which are continuously differ entiahle a.e. 
Then, with probability one, 

P(i?{„,}(Z(t)) I Z{t) = z) = hm P(i?|:]j(z) I Z{t) e A„(z)) 

can be calculated explicitly via the formulae given in Propositions \4-l\ \4-^ H-"^ '"^'^ 
RemarklJ^ 



Remark 4.3. We may also consider the case that G is countable. However, to transfer 
the results of the finite case, we have to ensure uniform convergence of the blurred 
intersection sets which is needed to compute V 1 Pf-Sl™'', (z)) in the 

denominator of Equation (|4.15p . To this end, we have to impose some additional 
conditions. For example, we could assume that for almost every z > there is only a 
finite number of shape functions involved in the intersection sets \A\ > 2. 

We are still left with simulating 112 | Z{t) — z given the occurence of a scenario 
£'{„^}(z) with J2a- ieA "-A = 1 for all i G {1, . . . , n}, that is, we are interested in 



fl {|n2n(CAx(o,oo)x{/})| = i} 



E{nA}i^) 



for Ga CR, f e G with {Ga x (0,oo) x {/}) n Ia{z) 7^ 0. Using Theorem with 
sets Ba^Ga-x (0,00) X {/}, A G 2{i'-"> \ we get that 

pf n {in2n/A(z) 5^1=1} £;{„.j(z)U hm n 



\A: nA = l 



A: nA = l 



Hi'a\^)) 



Thus, each random vector {Ta,Fa) G M x G, A e 2{i'- -"> \ 0, which is defined by 
|n2 n {{Ta) X (0, 00) X {Fa}) n lA{Z{t))\ = l if |n2 n lA{Z{t))\ = l, can be simulated 
independently. 

The distribution of {Ta, Fa) depends on the cardinal number of A. If j4 = {i} for 
some i e {1, . . . , n}, we have 

Pf({/}) 

P(TAeB,FA = /)-^^ ' y , B^B. 



24 



M. Oesting, M. Schlather 



For \A\ > 2, let /^(z) = {{tl^\ y^^\ h) , . . . , (^'-^ , /fe)}. Then, we get 

Pf ({/,}) 

-1 ter')=Ni/;'(t2-4''')-^2/;,(ti-t^^'')| 



Thus, we end up with the following procedure for calculating the conditional distri- 
bution of Z{tQ) given Z{t) with to,ti, . . . ,tn € z > 0. 

1. Compute the conditional probabilities (|4.15p for all the scenarios i?{„^}(z) and 
generate a random scenario following this distribution. 

2. For a given scenario £'{„^}(z) set 112 = {(T!a7 inin"^;^ -pji^^hrj) t Fa) '■ n-A — 1} 
where the law of (Ta, Fa) is given above. 

3. Independently, sample from lis = 11 n (R x (0, oo) x G) \ (-ft^(t,z) U J'^(t,z))- 
Then, Z{to) = max(^_„ j)gn2un3 "/(^o " «)■ 

In the next section, we will point out the capability of this exact approach by 
comparing it to other algorithms in the simple case of a deterministic shape function 
which is continuously differentiable. 



5. Comparison with the algorithm for the maix-linear model and 
transformation to Gaussian marginals 

Recently, Wang and Stoev [H] proposed an algorithm for exact and efficient condi- 
tional sampling for max-linear models 

Z{ti) — max aijYj, i = 0, . . . , n, 

where Yj, j — 1, . . . ,p, are independent standard Frechet random variables. 
Rewriting Z from p.ip as an extremal integral (see |20| ) 

Z{t) = Pf({/}) 7 fit - u)M^{du), t e K, 

/EG "^^ 

where Mi is a random sup-measure on R w.r.t. the Lebesgue measure, we can see that 
Z can be approximated arbitrarily well by a max-linear model, e.g. by 

ZnAt) = '^i^_^''\i_i^p^^-f3^^ ■ -^j yy + ^) ^) ■ ^ e h>o, 

where G = {/i, . . . , fk}- Then, we have ZM.h(t) Z(t) for any t € 'R as M — > cxo, 
/i ^ 0. 

We also consider another approach based on the assumption of a multi-gaussian 
model (cf. [2], p. 381). The data are transformed such that the marginal distribution is 
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Gaussian. As the marginals of Z are known to be standard Frechet, the corresponding 
transformation is given by 

: (0, oo) ^ R, $"^($i(a;)), 

where $ is the standard normal distribution function and $i — exp(— is the 
standard Frechet distribution function. The transformed random field Y — '^{Z) is 
stationary and second-order. As the covariance function C oiY can hardly be com- 
puted for general shape functions fi,. . . ,fk, we estimate it using maximum likelihood 
techniques, for instance, from a large parametric class like the Whittle-Matern class, 
i.e. 

assuming that y is a Gaussian random field. Under this assumption, the conditional 
distribution can be sampled easily (see [T3], for instance). Afterwards, the sample has 
to be retransformed via 

(0,00), y^<i>^\^y)). 

To compare these different methods, we need a measure for the goodness-of-fit of 
a distribution. Here, we use the continuous ranked probabilty score (CRPS) which is 
defined as [lO] 

/oo 
(Fi(j/)-l{^>,})2d2/, 
-00 

where Fi is a cumulative distribution function and a; € M. Note that CRPS{Fi,F2) := 
J CRPS{Fi,x)F2{dx) is a strictly proper scoring rule, i.e. 

CRPS{F2,F2) > CRPS{Fi,F2) 

for all cumulative distribution functions Fi, F2. If Fi and F2 both belong to measures 
with finite first moment, equality holds if and only if Fi = F2. Assuming that Fi has 
a finite first moment, the CRPS can be calculated via 

CRPSiFi,x) = l^F^ \X-X'\-Ef, \X~x\, (5.1) 

which shows CRPS{Fi,Fi) ^ -^Ep, \X-X'\ < 0. Here, X, X' - Fi are independent 
random variables. 

In order to compare different algorithms for getting a realisation from the conditional 
distribution 2'(fo) | Z{t), we consider K samples Zi, . . . , Zk of the random field Z. 
For each method m, we get an empirical distribution function F^™''^ as the (approx- 
imated) conditional distribution of log(Zi(to)) | ■^z(t), i = l,...,i4r, and calculate 
CRPS{Fl™'\\og{Zi(tQ))) via (j5.1|) . Here, we do the log-transformation to Gumbel 
marginals to ensure that the conditional distribution has finite expectation. 

Then, a measure for the goodness-of-fit is given by the mean score [lOj 

1 ^ 

CRPSK,m = - ^Ci?P5(^;^'"\log(Z,(to))). 
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Figure 3: Left: Construction of Z. The grey dots represent the points (s,ii ■ /(O)) with 
{s,u,f) £ n, the black hne is one reahsation of Z. The black dots mark Z{t). Right: 
Construction of Z conditional on Z{t). 



Gaussian extreme value process 



Conditional Gaussian extreme value process 





Furthermore, we have a look at the mean absolute error of the logarithmic conditional 
median 

1 ^ 

MAEK.,m = j^J2 ((^^"^'(0.5)) - log(Z,(io)) 

1=1 

Here, for computational reasons, we choose the Gaussian extreme value process [19) 
which has the deterministic shape function 



fix) = (p{x) = 



1 



exp 



1 . 

-r 



Furthermore, let n = 4, t = (—2,-1,1,2) and to — 0. Figure [3] shows two 
realisations of Z(-), the first one is sampled unconditionally and the second one is 
based on conditional sampling of the first one. 

The conditional distribution is calculated based on a sample of size 100 simulated in 
R [n]. The performance is measured via CRPSk,i / MAEka (conditional sampling 
via the Poisson point process), CRPSk,2 / MAEk,2 (conditional sampling for a max- 
linear model with M = 5 and h = 0.1 using the R package maxLinear [21]) and 
CRPSk,3 I MAEk.3 (conditional sampling via transformation to Gaussian marginals) 
with K — 1000 samples. As already mentioned, the latter approach requires the 
knowledge of the covariance structure of the transformed random field. This is assessed 
by first simulating data from this model on a dense grid repeatedly and then estimating 
the parameters of a Whittle-Matern covariance model based on maximum likelihood 
techniques implemented in the R package RajidomFields [18] . The parameters are 
chosen such that of the first and second method have a similar running time, the last 
method runs much faster. 

The results of the simulation study are shown in Table [TJ Here, CRPSk,i and 
MAEk,i can be interpreted as reference values as the first method is exact. We 
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Table 1: Results of the simulation study for f{x) = 75(2;) and K — 100. 





m = 1 


m — 2 


m — 3 


CRPSK.m 


-0.149 


-0.337 


-0.268 


MAEK,m 


0.217 


0.476 


0.357 



note that conditional sampling for max-linear models performs worse than conditional 
sampling via transformation to Gaussian marginals. 

For further analysis and comparison of these methods we do not restrict ourselves 
to pointwise prediction, but have a look at the sample paths. Additionally, pointwise 
quantile estimation of the conditional distribution can be done including the special 
case of the conditional median which can be seen as an analogon to Kriging. In case 
of conditional sampling via the Poisson point process and conditional sampling of a 
max-linear model the quantiles have to be estimated from the empirical conditional 
distribution. For sampling via Gaussian transformation the quantiles can be calculated 
directly by the means of Kriging result and variance. 

Figure 0] shows five sample paths and the median of the Gaussian extreme value 
process conditional on 

a. observations at four locations —2, —1, 1, 2, 

b. observations at eleven locations —2.5, —2, . . . , 2, 2.5. 

In general, conditional simulation via the Poisson point process yields sample paths 
which capture the main features of the process quite well. Even in the case of four obser- 
vations parts of the sample path are reconstructed exactly with a positive probability. 
For eleven observations most of the sample path is restored with high probability. 

The results of conditional sampling of the max-linear model are similar to the first 
method in case of four observations. For eleven observations, however, the method fails 
because of model misspecification. As the data to not match the discretized model, 
some observations cannot be reconstructed. For some realisations of the Gaussian 
extreme value process this problem even occurs in case of four observations. This is 
the main reason for the bad results of this method in the simulation study above. 

Gaussian transformation yields conditional sample paths which are structurally very 
different from the true ones. However, for eleven observations the deviations from the 
original sample path are quite small. 



6. Approximation in the case of an infinite number of shape functions 

Here, we drop the assumption that G is finite. Note that the measurability of II2 
and II3 is still an open question if G is uncountable. We present an approximation of 
the distribution of Z{tQ) given Z(t) based on a finite number of shape functions. 

Let Fi,F2, ■ ■ ■ be independent copies of F where F is defined as in Section [TJ Then, 
given Fi , . . . , F/v , we define 

ZN(t)= max uf(t-s), t€W^, (6.1) 
(s,«,/)en(«) 
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Figure 4: Comparison of the Gaussian extreme value process with different types of 
conditional simulations: a. simulations conditional on four observations at —2, —1, 1, 2, 
b. simulations conditional on eleven observations at —2.5, —2, . . . , 2, 2.5. In both cases the 
original Gaussian extreme value process (top left), conditional samples via the Poisson point 
process (top right) and conditional results for a max-linear approximation (bottom left) and 
an approximation by Gaussian transformation (bottom right) are shown. Black crosses: 
observations, grey lines: conditional sample paths, black line: conditional mean. 
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where H^^^ is a Poisson point process on R"^ x (0, oo) x {-Fi, . . . , Fn} with intensity 
measure 



1 



1 



A{AxBx{Fk}) = — I I —duds, A € B"^, B € B n {0,oo), k € {1 . . . ,n}. 



AJB 



Theorem 6.1. For any z > it holds 

n(^) I z{t) < z A n I z{t) < z 

as N ^ oo. In particular, ZN{to) | •^Ar(t) < z A- Z{to) \ Z{t) < z. 
Proof. We note that it suffices to show 

Urn P(|nW nMj\=nj, j = l,...,l\ Zjv(t) < z) 



p(|n n Mj\ = nj, j = i,...,i\ z{t) < z) 



(6.2) 



for z > and sets Mj = {aj,bj) x {cj, oo) x Gj c K"^ x (0, oo) x G which are pairwise 
disjoint. First, we consider the distributions of Z and Z^, respectively. It holds that 
P(Z(t) < z) = exp(-Eif), where 



H 



= / max 



F{ti - 1) 



dt. 



We define Hi,. . . , Hn replacing F hy Fi, . . . , F^, respectively. Then, we have 

P(^jv(t) <z) = E(P(Zw(t) <z\Fi,... .Fm)) = E l^exp ^iJ^^^ . 

As the sets Mi, . . . ,Mi are pairwise disjoint, the common distribution of 11 and Z{t) 
can be rewritten as 

p(|nnMj| = nj, j = i,...,i, z{t) < z) 
= p(Z(t)<z).I]p(^nnM,n|(x,2/,/)e5: y < , min^-^^^ij— ^| =n,^ 
= P{Z{t) < z) 



n n < (x, y, /) G {aj,bj) X [cj, oo) x Gj : y < min 



=l,...,n f{ti - x) 



Thus, we get 



p(|nnM,| = n,, z(t)<z) 

exp(- EH)Y[ ^ ^ / exp(- E H^'^), 



where 
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Analogously, we can define i/^.-'' replacing F by Fk in the definition of iJ^^^ for j — 
k^l,...,N. This yields 

P(|nWnAfj| =n„ Z{t)<z) 
- E („p (-1 £ H.) n -L (1 £ ^<...) "' „p (-1 g H'^') ) . 

AsHi,...,Hn ^^.^.d. H and h['\ iJ^^ for j = 1, . . . , Z with < 

'^"^i J- < oo and i?*^^-* < ^'=^''||^' — the Strong Law of Large Numbers gives 

jfJ2k=iHk EH and ^ Ef=i ^ Eij(j) a.s. as iV ^ oo. Thus, dominated 
convergence yields (|6.2I) and P(ZAr(t) < z) — P(Z(t) < z), which verifies the first 
assertion. 

The second assertion follows immediately by rewriting 

{Z{to) <z} = {|n n {{x, y,f)eS: yfito - x) > z}\ = 0}. 

If G is countable, we have 

lim nf ) I (ZAr(t) = z) = hm {Hn \ (i^t.z U i?^)) 

= n \ (i^t.z U i?^) = Ha I (Z(t) = z), 

where we used the first part of Theorem lG.ll and applied the second part of Theorem l3.2l 
to the processes Z and Zjy- This motivates to improve the approximation H « II^^^ 
by n ss Hj^' U Ha, i.e. by the following procedure: 

1. Simulate H^^^ | Z(t) = z. 

2. Independently of XIj^' , sample Ha | Z(t) = z which is defined as Iln (M x (0, oo) x 
G) \ (i^(t,z) U ^(t,z)) analogously to the second part of Theorem 13.21 

Then, Z{to) « max(5_„ j)gn2un3 uf{to - s). 

7. Application to the Brown-Resnick process 

We will apply the method of conditional sampling via the Poisson point process to 
the Brown-Resnick process |T| . 

Let {Wx{t), t e R}, X € (0,oo), be independent copies of a standard Brownian 
motion and — independently of the W^'s — let n be a Poisson point process on (0, oo) 
with intensity measure du. Then, 

Z{t) = max (^u exp (^W^{t) " y ^ ^ > t e 

defines a stationary max-stable process with standard Frechet margins. 

Recently, this process was generalized |12) and its mixed moving maxima represen- 
tation was given explicitly 0. 
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This is, 

Z{t)^ max M • iexp(-/(t - s)), i e M, (7.1) 
(s,M,/)gn 2 

where 11 is a Poisson point process on M x (0, oo) x C(M) with intensity measure 
ds du Pfl;(d/) and P^; is the law of the process 

R[t)^lt<oRi{-t) + lt>Mt). 

Here, i > 0}, {i?2(i), t > Q} are independent Bessel processes of a three- 

dimensional Brownian motion with drift ^ in its first component (cf. [151), i.e. 



Riit) ^ R2{t) ^ V(Wi{t) + |t|/2)2 + W2{t)^ + Ws{t)^, 

where Wi, W2 and W3 are independent standard Brownian motions. 

We will use the results obtained in the section above to sample from the conditional 
distribution of the Brown-Resnick process. However, the sample paths of exp(— _R(-)) 
do not satisfy the assumptions of Propositions 14.11 14.21 and 14.31 In particular, the 
sample paths are not continuously differentiable almost everywhere. To overcome this 
drawback, we do not use the exact sample paths F{-) ~ exp(— i?(-)), but the sample 
paths evaluated on a grid and interpolated linearly in between. 

We show that this procedure is correct in the limit. Let T — {t^, 2; G Z} C M with 
. . . < t_2 < t-i < to < ti < t2 < . . . such that limz_j._oo = —00, lim^-^oo — co 
and ||r|| = supj,g2(^z ^ tz-i)- Let {{t,FT{t)), t £ R} he the polygonal line through 
the points {{t,F{t)), t e T}. Furthermore, define Zrit) as in (|7.ip . replacing F by 
Fx • 



Proposition 7.1. For ||r|| — ^ we have Zxit) — > Z{t) in probability for all t e M. 
In particular, 

Zrito) I (^T(t) eB)^ Z{to) I {Z{t) e B) 
for all Borel sets B C M" with P{Z{t) e B) > and P(Z(t) e dB) = 0. 
Proof. W.l.o.g. we assume ||T|| < i. Then, we have 

f^FT{t)dt<J2 sup ^Frit) 

= sup Fiz + i) = exp ( - inf R{z + ) . 



, te[z-l/2,z + l/2] 

< > ; ( sup ^-F{t) + sup \F(t) 

te[z,z+i] ^ 
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We show that the rhs is L^-integrable. To this end, we assess 
E I Vexp (- inf R(z + t)] ] 

< 2 lim sup ^ E ( exp 

< 2 lim sup ^ E I exp 

fc=0 ^ 

< 2hnisup^E [ exp ( - W{k) + 



inf R(k + t) 

te[os] 



inf 

te[OA] 



k=0 



W{k + t) 
k 



k + t 



sup 

te[o,i] 



W{k + t)~W{k) + 



< 2 E exp sup 
V \telos] 



Wit) + - 



hm sup s 1 + E I exp 



k=l 



W{k) + - 



where is a standard Brownian motion. 
By [13] there exists some e > such that 



E \ exp I e sup 



W{t) 



< oo. 



Since \X\ < ^ +eX^, this impUes that E(exp(suptg[o l^i^) + |l)) is finite. 

Furthermore, using that exp (— ^a;^) da; < ^ exp (—-f-^ for aU xq > 0, we have 



E exp 



W{k) 



k 



2 

: exp 



-fe/2 
1 

/27rfc 



fc/2 



exp 



,2 

27rfc ' V 2fc 
(a; + kf 



-fc/2 ^ / ^ 

gj:+fe/2^^_ ( ) (Ij, 

V2^ V 2fc, 



2fc 



da; + exp(fc) 



-fc/2 



exp 



(x-k)' 
2k 



dx 



= , / exp dx + expl/e) / exp 

^/2^ V 2 y ^73^ 

2 / fc\ 2 / 9A: 

< exp — -J + exp(fc) exp 



V27rfc 
All in all, we get 



3v27rfc 



da; 



— . exp , 
3^2^ V 8 



E exp ^- ^ inf i?(z + j j < 2C lim sup + ^^^P 



< oo. 
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As we have an integrable majorizing random variable, dominated convergence yields 

dt 



lim E 

\\T\HO 

E I lim 

VllTIKO 



lFT{s-t)-^F{.s-t) 



lFT{s-t)-^F{s-t) 



dt] =0, s e K, (7.2) 



as F is continuous and integrable. 

By Theorem 2.1 in [20j Equation (L^ is equivalent to Zt(s) A Z{s) for all s G M. 

As 

m 

P(||(Zt(si), . . . , ZT{s.m)) - {Z{si), Z{s„,m >e)<J2 P(I^t(s.) - Z{s,)\ > e), 

i=l 

we get that all finite-dimensional marginal distributions of Zt converge in probability 
and therefore also weakly. Thus, for Borel sets A C M, B C R" with P{Z{t) e B) > 
and PiZ{to) e OA) = P(Z(t) e dB) = 0, we get 

lim FiZrito) G A \ Zrit) e B) = P(Z(to) G A \ Z{t) e B) 

||T||\0 

which completes the proof. 

Still, for any fixed T, the range of ^Ft is uncountable. Therefore, we have to use 
the approximation introduced in Section [6] 

We compare this approximation to conditional sampling based on the approach of 
[22] and the Gaussian marginals approach, see Section[SJ To compare these procedures, 
we simulate K — 100 independent samples of Z on the set {to, ^i, ^2, is, ^4} with = 0, 
ti = —2, <2 = —1, ^3 = 1 and = 2. We calculate the CRPS by sampling 50 
times from the (approximate) conditional distribution of log(Z(to)) given Z{t). Let 
CRPSK.ia I MAEkm and CRPSk^u / MAEk^u denote the CRPS / MAE based 
on two variants of conditional sampling of the Poisson point process for K samples, 
CRPSk.2 I MAEk,2 the CRPS /MAE of the approach by Wang and Stoev ^ 
and CRPSk.z I MAEk,3 the CRPS / MAE based on transformation to Gaussian 
marginals. For the Poisson point process approach, we chose N — 1000 as the number 
of shape functions on the grid T — {—7, —6.95, —6.9, . . . , 6.9, 6.95, 7}. However, if we 
restrict ourselves to a finite number of shape functions, the intersection set /b(z) with 
\B\ > 3 is most likely empty, even though Z{ti), i ^ B, may be determined by the 
same (s, u, f) G 11. Therefore, we do not only consider "exact" intersections, but also 
intersections which occur if the function values differ up to a given tolerance, i.e. we 
assume {t,y,f) G /b(z) with y = minj=i^..._„ 7(u^ 

< mm{y + £,y{l + e)} -4==> i^B. 



fiu - 1) 



for some given tolerance £ > 0. The simulation study is performed for e = 10 ^ 
yielding CRPSkm I MAEkm and e = lO'^ {CRPSK,ib / MAEK,ib)- Thus, by 
construction CRPSK,ia / MAE^^ia does not contain any intersections of more than 
two curves, but CRPSK,ib / MAEK,ib does. For the Wang-Stoev-approach, we chose 
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Table 2: Results of the simulation study for the Brown-Resnick process with A'' = 1000 and 
K = 100. 





m ^ la 


m = lb 


TO = 2 


TO = 3 


CRPSK,m 


-0.35 


-0.43 


-0.35 


-0.35 


MAEK,m 


0.49 


0.57 


0.49 


0.48 



M = 5, h = 0.05 and the same N — 1000 shape functions. The parameters are chosen 
such that the first and the second method have similar running times. 

Then, we get the results presented in Table [21 Note that the results are qprecise as K 
is small. However, for different runs we always got a similar behaviour. All the methods 
have almost the same accuracy, at least if e = 10~^. However, the last method runs 
much faster than the others. Furthermore, we notice the difference between CRPSkm 
and CRPSK,ib indicating that considering approximate intersections of at least three 
curves yields worse results. This is because these intersections involve incorrect shape 
functions. Furthermore, they lead to degenerated conditional distributions which are 
not supposed to occur in the case of the Brown-Resnick process. Thus, in this case the 
approximation of the mixed moving maxima process seems to be appropriate only if 
we do not consider intersections of three or more curves. 



8. The discretized case 

By now, we have considered the general model ()l.ip . The procedure we proposed is 
exact in the case of a finite number of shape functions which are sufficiently smooth. 
However, as the example of the Brown-Resnick process in Section [7] illustrates, we may 
run into problems if these assumptions are violated. 

Now, we modify our general model (II. ip and use a discretized version 

Z(t) = max uf{t-s), tepZ'^, (8.1) 
(s/it,/)gn 

where H is a Poisson point process on pZ"^ x (0, oo) x G where p > and G C (0, oo)p^ 
is countable. The intensity measure of H is given by 

A{{s} xBx {g}) = J2 -^P^ds) X / u'^ du x PpUg}) 

where Pi? is a probability measure belonging to a G- valued random variable F with 
E(E.ez- F{pz)) = 1. 

Using the same notations as before, we get the same results as in Section[3l However, 
all the calculations can be done explicitly without any further assumptions on / e G. 
We get the following results. 

Proposition 8.1. Let B ^ {i} C {1, . . . ,n} , z £ and 

Di{z) = {{xj) epZ'^ xG: ix,yj) G /{,}(z) for some y e M}. 
Then, we have 

{xj)eDi{z) 
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Proposition 8.2. Let B e 2{i'- -"> \ 0, \B\ > 1 and z e K" such that 

IaIz) ^ {{xj,yj,f.j), j = l,...,/(z)} 
with l{z) > 0. Then, for m large enough, we have 

In particular, A(4"^(z)) G 0(2-™), &ut A(4'"^(z)) ^ 0(2-^(1+'^)) for any e > 0. 

By these formulae, all the scenario probabilities can be calculated. As the intensity 
of each intersection set has the same rate of convergence, only scenarios with minimal 
|n n Kt^z\ occur. 

We note that our model is very close to the model investigated by [22]. To see this, 
we calculate that 

P[Z(ti) < zi, . . . ,Z{t„) < Za) ^ exp \ - } } max 

\ ^-^ i=l,...,n Zi 

Therefore, we get that 

Z ^ max max (/(• - pz) Ppi{f})z'fA , 

where the random variables z <^U^, f £ G, are independently standard Frechet 

distributed. 

This means, the model (j8.1l) is a max- linear model and the exact conditional distri- 
bution can be calculated via the algorithm of [52] if G is finite and the support of each 
/ G G is finite. However, the algorithm of conditional sampling of the Poisson point 
process and the algorithm of [22] do not work in exactly the same way. According 
to the latter algorithm one samples from each random variable Zj. . This procedure 
corresponds to simulating the largest point of 11 n (0, oo) x {pz} x {/} for each z g Z, 
f G G. The first algorithm includes the simulation of points in 11 until a terminating 
condition given in Theorem 4 of [17) is met. Computational experiments show that 
both algorithms yield identical results. Also the technical results are related. 

For example, the occurence of a scenario J C Z'' x G (in the notation of Wang/Stoev) 

z(t - ^ 

corresponds to the event that 112 consists of |J| elements {pz, mini^i^...,„ f(ti-pz) ' 
with {z, f) G J. By this correspondence, the statements 

• |n n Kt^zl is minimal a.s. 

• an occuring hitting scenario J satisfies |J| — r{J'{A,x)) a.s. ([22]) 

are equivalent, both claiming that the number of points generating the observation 
(t, z) is minimal. Hence, although the approaches look quite different, there are similar 
observations and results in [22] and in this section. 
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